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The thermodynamic behavior of the two-flavor(A'^j^ =2) three-color {Nc = 3) Polyakov- 
loop-extended Nambu-Jona-Lasinio model at the finite chemical potential is investigated. 
New lattice gluon data for gluon thermodynamics are used defining the effective potential 
within polynomial and logarithmic forms of its approximation. We study the effects 
of using different sets of data and different forms of the potential on thermodynamic 
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attention is paid to the phase diagram in the (T, fi) plane. 
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1. Introduction 

In the last years the Polyakov-loop-extended Nambu-Jona-Lasinio (PNJL) 
model was widely used in the study of thermodynamics and the phase diagram 
of hot and dense matter. Results of this research are expected to play an impor- 
tant role in our understanding of the evolution of the early universe and physics of 
heavy ion collisions at relativistic energies. This improved field theoretical model is 
fundamental for interpreting the lattice QCD data and extrapolating into regions 
not yet accessible for lattice simulations. 
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An attractive property of the PNJL model is the synthesis of the Polyakov 
loop dynamics with the Nambu-Jona-Lasinio model, combining the two principal 
nonperturbative features of low-energy QCD: confinement and spontaneous chiral 
symmetry breaking. A particular feature of this model is that one can uniquely 
determine the coupling between the chiral condensate, which is an order parameter 
of the chiral phase transition when -^0, and the Polyakov loop, which is the 
order parameter for the deconfinement phase transition in the limit rriq — > oo. The 
model is re markably successful in reproducing lattice data on the QCD thermody- 
namics Etzmni 

However, the choice of the parameter set as well as regularization of integrals, 
as noted by several authors, is a nontrivial question. As it is well known, the order 
of the phase transition in the (T, /z) plane is sensitive to the parameter choice. It 
was already noted SEHl^jiat one can choose different sets of parameters which allow 
for a first order phase transition, giving a reasonable fit to physics observables in 
the vacuum but predicting different physical scenarios at finite temperature T and 
chemical potential fi. In addition, recently new lattice data for the pure gluon QCD 
sector defining the effective Polyakov loop potential have been obtained which 
differ noticeably from the old data^S. 

In this paper we investigate how the input information from lattice QCD and the 
used forms of the effective potential influence general properties of thermodynamics 
at finite temperature T and baryon chemical potential fi. After introduction, in Sect. 
2 we consider the polynomial and logarithmic parameterizations of the Polyakov 
loop effective potential for the new and old pure gluon lattice data within the two- 
flavor PNJL model. Independent of the temperature, the model parameters defined 
by properties of quarks and mesons are presented in Sect. 3. Comparative study of 
the thermodynamics and phase structure, their dependence on the lattice input and 
used parametrization are considered in Sect. 4 at finite T and /i. The last Section 
summarizes the obtained results. 



1.1. The Nambu Jona-Lasinio model with Polyakov-loop 

The deconfinement in the pure SU{Nc) gauge theory can be simulated by introduc- 
ing an effective potential for a complex Polyakov loop field. The PNJL Lagrangian 
employed in this work is 

/:pNJL = g(^7M^"-™o)<Z + G[(gg)' + (g^75r<z)'] -W($[A],$[A];T) . (1) 

Here, a local chirally symmetric scalar-pseudoscalar four-point interaction of quark 
fields q, q is introduced with an effective coupling strength G, f is the vector of 
Pauli matrices in fiavor space, toq is the diagonal matrix of the 2-flavor current 
quark masses, mp = diag (m^, m^), = to^ = ™o- 

The quark fields are coupled to the gauge field A'^ through the covariant deriva- 
tive Z?^ = 9^ — iA^. The gauge coupling g is conveniently absorbed in the definition 
A^'{x) = gA'^,^ where A'^ is the SU{'i) gauge field and Aa is the Cell-Mann ma- 
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trices. The gauge field is taken in the Polyakov gauge A'' — 5qA^ ~ —iS'^Ai. The 
field $ is determined by the trace of the Polyakov loop L{x) and its conjugate^ 

= ^Tr,L(f) , m = ^Tr,Lt(f) , 



where L[x) — Vexp 



dTAi{x, t) 



(5 — 1 /T being the inverse temperature. In 



the absence of quarks, we have $ = $ and the Polyakov loop servers as an order 
parameter for deconfinement. 

The gauge sector of the Lagrangian density ([T]) is described by an effective 
potential W T) . The effective potential must satisfy the Z{3) center 

symmetry. In accordance with the underlying Z{3) symmetry, one can choose the 
following general polynomial form: 

y4 ^ =-^$$--| ($^ + $^)+^ (2) 



b,{T)=ao + a,^^j+a,^-^j + as ^-^ j . (3) 

The Z{3) symmetry leads to some freedom in the choice of the effective potential 
form. Along with the simplest polynomial form, Eq.Q, there exists an expression 
with a logarithm in place of the higher order polynomial terms in $ 1^. In the 
logarithmic form the potential is 

^ ' = --a (T) + 5 (T) In [l - 64>$ + 4($3 + $3) _ 3($$)2] ^ (4) 

a (T) = ao + 5i (^^) + 52 ' , 6 (T) = 63 ) ' . (5) 

The pressure of a pure-gauge system is given by p = —U. 

Pure-gluon lattice data from Ref. are tradit ionally used to find the parameter 
set for both forms of the effective potential I 3 | 23 | jj^ contrast, our work is based on 
new gluon lattice data looking for a new potential to fit the lattice pressure. In 
finding the potential parameters the following conditions should be satisfied: $ -> 1 
and p/T'^ — > 1.75, when T — >■ 00. As immediately follows from these conditions, 
ao = 3.51 for the logarithmic potential and the constraint 1.75 = ao/2 + 63/3 — 64/4 
for the polynomial potential. Minimizing W($,$,T) with respect to variation of 
$ and taking into account that $ = $ at = 0, we can find parameters using 
the method of least mean squared deviations. Thus, for the critical temperature 
To =270 MeV the following parameter sets were obtained (see Table [T] and Table 

In Fig.[Tl old and new lattice gluon data are compared together with the results of 
their approximations. As is seen the lattice results differ by about 10% at T/Tq > 2. 
New data are plotted by circles but the density of measured points is so high that 
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the results look like a shaded band. Both polynomial and logarithmic forms are in 
nice agreement with the data and it is hard to distinguish them from each other. 

In general, the^arameter Tq depends on the number of active fl avo rs and the 
chemical potential ^ . In the pure gauge sector Tq = 0.27 GeV was use dEI. The effec- 
tive potential for both sets of parameters at T = 0.2, 0.27, 0.32, 0.54 GeV is shown 
in Fig. [2l Both sets describe quite satisfactorily the Polyakov loop as a function of 



Table 1. Parameters of the effective potential U[A] with 
the polynomial form. 



ao ai 


112 


as 




b4 


old dataESi 6.75 -1.95 


2.625 


-7.44 


0.75 


7.5 


new datal^ 6.47 -4.62 


7.95 


-9.09 


1.03 


7.32 


Table 2. Parameters of the effective poten- 




tial U "I>; r) in the logarithmic form. 








ai 


5,2 






old dataEU 3.51 


-2.47 


15.2 


-1.75 




new datal^ 3.51 


-5.121 


20.99 


-2.09 





1.0 




Fig. 1. Scaled pressure in the pure gauge sector as function of scaled temperature. The old 
and new^^ lattice data are plotted by circles and triangles, respectively. Solid lines correspond to 
polynomial form of potential and dashed lines correspond to the logarithm form. 
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Fig. 2. The Polyakov loop effective potential W as a function of <& for various values of temperature 
for old (top) and new (bottom) sets of parameters. Left panel corresponds to the polynomial form 
and right panel corresponds to the logarithm form of the potential. 



temperature. In accordance with the Z(3) center symmetry, the following properties 
of the effective potential U{^, $; T) are seen. At low temperature W(<I>, T) has a 
single minimum at $ = (a confinement phase); the effective potential is getting flat 
for the critical temperature T = Tq and above critical temperature (a deconfinement 
phase) a second minimum arises at nonzero $, as a consequence of Z(3) symmetry 
breaking; in the T ~> oo limit, $ — s> 1 (see Fig. One should note that after the 
second minimum the logarithmic potential $ increases faster than the polynomial 
one forming a more distinct minimum. With the introduction of quarks the critical 
temperature goes down. The range of applicability of this model is T < 2.5Tc since 
at higher temperature transverse gluons start to contribute significantly 

The Polyakov loop is compared with the lattice results in Fig. [3l In reasonable 
agreement of both forms with the lattice data, the new parameter set predicts 
slightly lower values of the field $ because the pressure for new data is also below 
the old one. 



2. Quarks and light mesons in the PNJL model 

The grand potential density for the PNJL (Nt =2) model in the mean-field approx- 
imation is given by the following equation '^^'^ 
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Fig. 3. Temperature dependence of the Polyakov loop $ for polynomial ^ (left) and logarithmic 
l l4t (right) forms. Lattice data are from Solid and dashed lines correspond to new and old 
parameter sets, respectively. 



where the quark term is 



(6) 



-2N,Nf 



and the functions are 



E„ - 2NfT 



J [111 ^* (^P) + In ^* (^f )] 



1 + 3 $ + '^e^P'^p e 



(Ep) 1 + 3 ( $ + $e~^^p ) e-^^p + e 



(8) 
(9) 



where Ep — y^p^ + mn? is the quasiparticle energy of the quark; E^ — Ep ^ the 
upper sign applying for fermions and the lower sign for antiparticles. 

Since NJL-type models are nonrenormalizable, it is necessary to introduce a 
regularization, e.g., by a cutoff A in the momentum integration. Following we 
use in this study the three-dimensional momentum cutoff A for vacuum terms and 
extend this integration to infinity for the matter contributions given by the second 
term of Eq. [71 A comprehensive study of the differences between the two regulariza- 
tion procedures (with and without cutoff on the quark momentum states at finite 
temperature) was performed in 

In the mean-field approximation, we can obtain the constituent quark mass m 
from the condition that the thermodynamic potential ([6]) will have a minimum with 
respect to variati on of this parameter, dfl/dm = 0. This condition is equivalent to 
the gap equation I 9 I ^6 | 



m = mo — 2G (qq) , 



(10) 
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where the quark condensate is defined as {qq) = dfl/drriQ. For the mass gap equation 
we get 



m 



m = mo+4GN.N,J^^-[l-f--r] (11) 
with the modified Fermi-Dirac distribution functions for fermions and antifermions 

^-PEl j_ ,-3/3E+l 



/+= [^ + 2^6-"^^ )e'P^p +e-^^^p /N^iEp) , (12) 



r 



!> + 2$e-'^^p") e-^'^p + e-3^^p"] /N^{Ep) . (13) 



Moreover, for PNJL calculations we should find the values of $ and <& by minimizing 
fl with respect to $ and $ ^ at given T and /i. One should note that if $ 1, the 
expressions Eqs. p^ . p^ reduce to the standard NJL model. 

For a self-consistent description of the particle spectrum in the mean-field ap- 
proximation, the meson correlations have to be taken into consideration. These 
correlations are related to the polarization operator of constituent fields. For scalar 

and pseudosca lar particles the polarization operators are represented by loop inte- 
gralsttOmiZl 

^afiP') - / Tr Hr'-Sip + P)^J,r'^S{p)] , (14) 

nff (P^) = I ^Tr [Sip + P)S{p)] , (15) 

where S{p) is the quark propagator and the operation Tr is taken over Dirac, fiavor 
and color indices of quark fields. 

From the point of view of the polarization operators, the pseudoscalar (tt) and 
scalar (a) meson masses can be defined by the condition that for = (M^a) the 
corresponding polarization operator 11'^^ (Af^) (H'^'^rM^)) leads to a bound state 
pole in the corresponding meson correlation function 121. For mesons at rest (P = 0) 
in the medium, these conditions correspond to the equations 

l + ^CVciV^/lfj^^g^d-r-D^O, (16) 

In order to solve Eqs. pU|) . (fTB|) and (fT7|) . a set of model parameters has to be 
determined: the above-mentioned cutoff parameter A, the current quark mass mo 
(in the chiral limit toq — 0) s-'id the coupling constant G. These parameters are 
fixed at T = to reproduce physical quantities: the pion mass Mtt = 0.139 GeV, 
the pion decay constant F^^ = 0.092 GeV and the quark condensate {qq)^^^ = —250 
■'^sV- The ob tained parameters are the same as in those obtained earlier in our 
papers I '^^ l '^^ l and are shown in Table [3l Meson masses obtained as solutions of the 
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Table 3. The set of model parameters reproducing ob- 
servable quantities (in brackets) and the chiral condensate 
(qg)i/3 = _250 MeV. 



mo [MeV] 


A [GeV] 


G [GeV] -2 


Fvr JGeV] [GeV] 


5.5 


0.639 


5.227 


(0.092) (0.139) 



gap-equation (|T0|) and Eqs. (fT6|) . (fTTl) at nonzero T are presented in Fig. |4] for two 
parameter sets. 




01 ■ ' ■ ' ^ — '-^=1 0.01 ■ ' ■ ' <- 

0.0 0.1 0.2 0.3 0.0 0.1 0.2 0.3 



T (GeV) T (GeV) 

Fig. 4. Temperature dependence of the masses rrig, and Ma masses at fi — GeV for 
polynomial (left panel) and logarithmic (right panel) forms of potential. The PNJL results for new 
and old parameter sets are given by the solid and dashed lines, respectively. The Mott temperatures 
for both parameter sets are plotted by vertical lines. 

The temperature modification of the quasiparticle properties is clearly seen in 
this figure. Up to the Mott temperature Tjviott, defined as M7r(TMott) = 2mg(TMott), 
the a mass practically follows the behavior of 2mq{T) with a drop toward the pion 
mass signaling partial chiral symmetry restoration. At T > Tiviott the masses of 
chiral partners become equal to each other, Afo- w M-j^. and then both masses 
monotonically increase with temperature. Below the Mott temperature, the pion 
mass remains practically constant. It justifies that TMott is a little bit lower for the 
new parameterization than for the old one when the polynomial form is used, while 
both values of TMott coincide in the case of logarithmic parametrization. 

3. Thermodynamics of the PNJL models 

The thermodynamics of particles is described in terms of the grand canonical en- 
semble which is related to the Hamiltonian H as follows: 



(18) 
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where N is the particle number operator, n is the quark chemical potential and the 
operator Tr is taken over momenta as well as color, flavor and Dirac indices. If O is 
known, the basic thermodynamic quantities - the pressure p, the energy density e, 
the entropy density s, the density of quark number n and the specific heat Cy - can 
be defined as follows: 



P = -T7 



e 
n 



df), 



p + Ts + iin, 



_T(ds\ 
V VST A 



(19) 

(20) 
(21) 
(22) 

(23) 



The thermodynamic potential in equilibrium corresponds to a global minimum 



1.0 

0.8 

0.6 

0.4 

0.2 

0.0 
1.0 

0.8 

0.6 

0.4 

0.2 



0.0 L— 
0.5 



polynomial form 




logarithm form 




T/T 




Fig. 5. The temperature dependence of the scaled pressure and energy density within the PNJL 
model for /i = in two schemes of parametrization. Dotted lines correspond to the old data 
parameterization and solid lines are new ones. Lattice data points for A^^ = 2 at /x =0 are from 
Ref. 30. Circles, squares and diamonds correspond to calculations at Nt =6 with the mass ratio of 
the pseudoscalar to vector meson mps/mv =0.65, 0.70 and 0.75, respectively. 
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with respect to variations of the order parameter 

dm ' dm'^ 

All these relations (HH) , ([201) , (l2l|) , ([22]) , ([231) describe the thermodynamics of the 
system. For the considered models the thermodynamic potentials are defined from 
Eq. ([6]) . From this equation we can read off the vacuum part 

n^ac = ^^"^ - 2N^Nf J ^E,. (25) 

This quantity does not vanish as T — > and /i — > 0. Therefore, in order to obtain the 
physical thermodynamical potential which corresponds to vanishing pressure and 
energy density at (T, /i) = (0, 0), one has to renormalize the thermodynamic poten- 
tial by subtracting its vacuum expression p5|) . This corresponds to the following 
definition of the physical pressure 

P p{T,fJ.,m) -p{0,0,m) 



ji4 2^4 



(26) 



With increasing temperature the pressure has to reach the Stefan-Boltzmann 
limit which in the chiral limit for the PNJL model is given as 

^ = (iV2 N^Nf— ~ 4.053 , (27) 

r4 ^ M5 M80 ' ^ ' 

where the first and second terms correspond to gluons and quarks, respectively. 

If the regularization A = 0.639 is used, the T-behaviour of the thermodynamic 
quantities considered is roughly the same while t heir absolute values are noticeably 
lower, being far from the Stefan-Boltzmann limit ^i^. 




Fig. 6. Comparison of the scaled quark number density as a function of temperature at fi =0.8 
and 0.6 GeV. Lattice data points are from Ref. 31. 
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Within the PNJL model with A — )■ oo (we can use A — >■ oo because most of 
the integrals in the PNJL are convergent) the reduced pressure and energy density 
exhibit reasonable behavior consistent with the recent lattice QCD results for the 
vanishing chemical potential ^221 (gee Fig.[S]) keeping in mind that the mps/mv ratio 
in lattice calculations is still far from that for physical masses mps/mv 0.2. As 
another example of thermodynamic characteristics, the temperature dependence of 
the reduced quark number density riq/T^ (see Eq. (|22p l is presented in Fig.[6l Model 
results are in good agreement with the lattice data for both values of the chemical 
potential considered. The logarithmic approximation of the effective potential U 
seems to describe lattice data better than the polynomial one. 

3.1. Phase diagram and the CEP 

Within NJL-like models there are several characteristic temperatures. The param- 
eter Tq entering into the effective potential ([3]) of the PNJL model has been noted 
above. For 7r-mesons, the Mott temperature TMott is provided by the condition 
MTr(TMott) = 2771, (TMott)- Above TMott the pion dissociates into a quark and an- 
tiquark and does not exist as a bound state. Similarly the a meson dissociation 
temperature Tj" is given by the equation M^(T^) = 2M^{T^) EES], other char- 
acteristics of phase transitions are the pseudo-critical temperature for the chiral 
crossover T-^, defined by the maximum of d{qq) / dT , and the pseudo-critical tem- 
perature for the crossover deconfinement transition Tp that can be found from the 
maximum of d^/dT. Their difference is less than 0.013 but it increases with de- 
creasing Tq. The third quantity is Tc assumed to equal T^. But in the PNJL model 
it is higher than the lattice result Tc ~ 192 GeV. Thus, it was suggested to define 
Tc as an average of two transition temperatures T^ and Tp ^^El 
All these quantities obtained a.t ft — are presented in Table 21 



Table 4. Characteristic temperatures in the PNJL models for fi =0. 







Tx 


Tp 


Tc 


^Mott 




Polynomial form of potential 


new 


0.2455 


0.2335 


0.2395 


0.259 


0.247 




old 


0.2575 


0.2485 


0.253 


0.27 


0.257 


Logarithmic form of potential 


new 


0.2305 


0.2295 


0.23 


0.264 


0.2645 




old 


0.2345 


0.2335 


0.234 


0.2645 


0.252 



To define the crossover transition line, the chiral condensate < qq > and the 

Polyakov loop $ were used as the order parameters. As shown in Fig. [71 these 

quantities are the temperature-dependent functions and demonstrate a quick change 

near the transition line which essentially depends on the chemical potential /i. The 

position of this line is defined by local maximum oid < qq > / dT =0 and d^ / dT =0. 

To find the first order transition line, it is convenient to introduce the baryon number 
dn„ 

susceptibility Xq = —T^\T=const- The first order phase transition ends just at point 
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where Xq has a pronounced maximum and this point is called critical endpoint 
(CEP) where the phase transition of the second order. At T > Tcep the baryon 
number susceptibility has a sharp rise and it can be considered as the presence of 
an ideal gas of weakly interacting quarks. 




T (GeV) T (GeV) 



Fig. 7. Temperature dependence of the chiral condensate and Polyakov loop potential with the 
new set of parameters within the PNJL model. Solid, dashed and dot-dashed lines are calculated 
for fi =0, n = t^CEP and fi > fiCEP, respectively. 

The behavior of the baryon number susceptibility Xq as a function of the chemical 
potential for three different temperatures around the CEP is presented in FiglS] 
For T < Tc EP , fJ- > fJ-c ep we have a phase transition of the first order with clear 
discontinuity; for T ~ Tcep the susceptibility Xg diverges at /i = hcep', for T > 
Tcep the discontinuity at the transition line disappears and we observe crossover 
type of the phase transition. The polynomial and logarithmic approximations of 
the Polyakov loop predict very similar results. Similar behavior exhibits also the 
specific heat c„. 

As is seen from Fig. |8l both models show the CEP at the temperature 
Tcep — T^ below which the chiral phase transition is of the first order. At 
this poin t (Tce p , Uc e p ) the phase transition changes from the first order to 
crossover 1 1 4 | 1 5 | 1 6 | ^^ns point the second order transition is present. 

The phase diagram in the (T, /i) plane is presented in Fig. |9l Within the PNJL 
model the positions of the critical endpoints (Tcep, fJ-CEp) are (0.118,0.3166), 
(0.11,0.3192) for the logarithmic form and (0.10,0.3175), (0.09,0.322) for the poly- 
nomial form, where the first pair of numbers correspond to the new d ata set and 
the second one is for the old data set (in GeV). As was noted in Ref. ^^^l^ critical 
properties of observables are significantly influenced by the chosen parameter set 
and regularization procedure. As follows from Fig. [9l the substitution of the new 
basic lattice data with using the polynomial and the logarithmic forms for the U 
approximation influences more significantly the chemical potential of the critical 
endpoint ficEP rather than its temperature Tcep- 
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0.1 0.2 

H (GeV) 



0.3 



Fig. 9. Phase diagrams of the PNJL model with polynomial (left panel) and logarithmic (right 
panel) forms of the potential. Solid lines correspond to the first order phase transition, dotted and 
dot-dashed lines are crossover for new and old lattice data, respectively. The phase region near the 
critical endpoints is zoomed in the insert. 



4. Summary and conclusions 

We have considered the PNJL {Nc = 3, Nf =2) model and investigated its phase 
structure at finite T and fj,. Two different sets of parameters based on the old and 
new lattice data for the pure gluon sector were used and for each of these sets the 
two different parameterizations of the effective potential ^/($,$,r) - polynomial 
and logarithmic - were applied. The thermodynamics for all the developed versions 
of the PNJL model was studied and compared with the available lattice data. Con- 
sideration of different thermodynamic observablcs like pressure and energy density, 
their T and /U behavior as well as the quark number density serves as an important 
probe of the model. We found that in spite of a noticeable disagreement between 
the old and new original lattice data, the effective gluon potentials U are quite close 
to each other and a larger difference is due to the form of their approximation: the 
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logarithmic form predicts a more distinct and narrower minimum at high T . 

The model qualitatively reproduces both tt and a meson properties in hot, dense 
quark matter and the rich and complicated phase structure of this medium providing 
information on the order of phase transitions and the position of critical points. This 
information depends stronger on the form of the effective potential rather than on 
the used lattice data set. 

Unfortunately, the position of the calculated CEP in the (T, ^) plane is still far 
from the predictions of lattice QCD and empirical analysis. Further elaboration of 
the presented model is needed. In particular, the inclusion of entanglement inter- 
actions between quark and gauge deg rees of freedom in addition to the covariant 
derivative in the original PNJL model and the incorporation of explicit diquark 
degrees of freedom are of great interest. Both modifications ^^^-'^ reproduce lat- 
tice data at /i > better than the original PNJL model and influence the position 
and the nature of the critical endpoint in the (T, /i) phase diagram. Moreover, with 
the use of the logarithmic form of the effective potential, these models result in the 
appearance of new phases. 

One should note that we are restricted to the case without diquark correlations 
and thus possible color superconducting phases at low T and high /i are ignored. 
It is attractive also to include into consideration the color superconducting phases 
and nonlocality of the interaction 1^31 as well as efiiects beyond the meanfield^^^^^. 
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